      subroutine read_para
      include 'cst.inc'
      
c	the following four files are obtained from the auxillary 
c	program setgrid
c
	open(unit=1,file='../SetGrid_FSU/gridpara.dat')
	read(1,*) imax,jmax,htop,hbot,aa,const,deltax,deltay
	close(unit=1)
c
	open(unit=1,file='../SetGrid_FSU/gshore.dat')
	read(1,'(i3,f7.2)')(jshore(i),yoff(i),i=1,imax)
	close(unit=1)
c
	open(unit=1,file='../SetGrid_FSU/shoreln.dat')
	read(1,*)
	read(1,*)
	read(1,'(3(f9.2,2x))')(xshore(i),yshore(i),sangle(i),i=1,imax)
	close(unit=1)
c
	open(unit=1,file='../SetGrid_FSU/dgrid.dat')
	read(1,*)
	read(1,*)
	do j=1,jmax
	   do i=1,imax
	      read(1,*)x,y,h(i,j),a(i,j)
	   enddo
	enddo      
	close(unit=1)

c
c	read in from the program input file para.dat
c
	open(unit=1,file='para.dat')
	read(1,*)
	read(1,*)acent
	read(1,*)
	read(1,*)astorm
	read(1,*)
	read(1,*)beta
	read(1,*)
	read(1,*)coef
	read(1,*)
	read(1,*)dcyyy
	read(1,*)
	read(1,*)dczero
	read(1,*)
	read(1,*)dcmax
	read(1,*)
	read(1,*)dlyyy
	read(1,*)
	read(1,*)dlzero
	read(1,*)
	read(1,*)dlmax
	read(1,*)
	read(1,*)dt
	read(1,*)
	read(1,*)edfact
	read(1,*)
	read(1,*)gfact
c	read(1,*)
c	read(1,*)hwave_max
c	read(1,*)
c	read(1,*)hwave_min	
      hwave_max=-30.0
      hwave_min=-3.0
	read(1,*)
	read(1,*)itmax
	read(1,*)
	read(1,*)k1
	read(1,*)
	read(1,*)k2
	read(1,*) 
      read(1,*)mfact
	read(1,*)
	read(1,*)nfact
	read(1,*)
      read(1,*)nnsave
	read(1,*)
	read(1,*)pore
	read(1,*)
	read(1,*)scr
	read(1,*)
	read(1,*)rfact
	read(1,*)
	read(1,*)wo
	read(1,*)
	read(1,*)qriv
	read(1,*)
	read(1,*)iriv
	close(unit=1)
c
	acent = ((4.*atan(1.0))/180)*acent
	astorm = ((4.*atan(1.0))/180)*astorm
	
	return
	end